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ABSTRACT 



Context. Recent works with improved model atmospheres, line formation, atomic and molecular data, and detailed treatment of 
blends, have resulted in a significant downward revision of the solar oxygen abundance. 

Aims. Considering the importance of the Sun as an astrophysical standard and the current conflict of standard solar models using the 
new solar abundances with helioseismological observations we have performed a new study of the solar oxygen abundance based on 
the forbidden [Oi] line at 5577.34 A, not previously considered. 

Methods. High-resolution (R > 500 000), high signal-to-noise (S/N > 1000) solar spectra of the [Oi] 5577.34 A line have been 
analyzed employing both three-dimensional (3D) and a variety of ID (spatially and temporally averaged 3D, Holweger & Miiller, 
MARCS and Kurucz models with and without convective overshooting) model atmospheres. 

Results. The oxygen abundance obtained from the [Oi] 5577.3 A forbidden line is almost insensitive to the input model atmosphere 
and has a mean value of log 6o = 8.7 1 ± 0.02 (cr from using the different model atmospheres). The total error (0.07 dex) is dominated 
by uncertainties in the log gf value (0.03 dex), apparent line variation (0.04 dex) and uncertainties in the continuum and line positions 
(0.05 dex). 

Conclusions. The here derived oxygen abundance is close to the 3D-based estimates from the two other [Oi] lines at 6300 and 
6363 A, the permitted O i lines and vibrational and rotational OH transitions in the infrared. Our study thus supports a low solar 
oxygen abundance (log eo ~ 8.7), independent of the adopted model atmosphere. 
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1. Introduction 

In recent years a significant reduction in the solar oxygen abun- 
dance has been proposed from log eo = ^og(No/Nu) + 12 = 8.93 
(Lambert 1978; Sauval et al. 1984; Grevesse et al. 1984; Anders 
& Grevesse 1989) down to log eo ~ 8.7 or even lower (Allende 
Prieto et al. 2001, 2004; Holweger 2001; Asplund et al. 2004; 
Melendez 2004; Socas-Navarro & Norton 2007; CafFau et al. 
2008; but see Ayres et al. 2006, Ayres 2008, and Centeno & 
Socas-Navarro 2008 for a departing conclusion). This large revi- 
sion in the oxygen abundance has an important impact in many 
areas of astrophysics. In particular, solar models computed with 
the low oxygen abundance do not agree with precise helioseis- 
mological measurements (see Basu & Antia 2008 and references 
therein). 

If the revised solar oxygen abundance is correct, this would 
imply that the standard solar structure model may be missing (or 
have a too simplified treatment of) important physical processes, 
which would carry over to errors in the calculations of stellar 
evolution at large. On the other hand, if the fault lies with the 
new solar photospheric abundance analyses it raises concerns 
about how well we understand stellar atmospheres and spectral 
line formation in general with wider implications for Galactic 
chemical evolution studies relying on accurate stellar abundance 
measurements. Indeed, the oxygen content in metal-poor stars 
is currently under debate (see Asplund & Garcia Perez 2001; 
Melendez et al. 2001, 2006, and references therein). 



Since the Sun is used as a fundamental reference in most ar- 
eas of astrophysics, it is important to study all spectral features 
available for abundance analysis. Recent studies have mainly 
iised the [Oi] 6300, 6363 A forbidden lines, the Oi 7777 
A triplet, pure rotation OH lines as well as the fundamental 
vibration-rotation OH lines (Allende Prieto et al. 2001, 2004; 
Asplund et al. 2004), from which log eo - 8.66 + 0.05 have been 
recommended (Asplund et al. 2004). Melendez (2004) presented 
a study of the first-overtone OH lines, which are unfortunately 
very weak for a precise determination of the oxygen abundance, 
but within the uncertainties these lines also support a low solar 
oxygen abundance. 



Another independent way to estimate the solar oxygen abun- 
dance is using the forbidden line at 5577.34 A. This line has 
been disregarded in the past due to blending with C2 lines (e.g. 
Altrock 1968; Lambert 1978). Here, we present the first detailed 
study of the [O i] 5577 A line, carefully taking into account the 
blends by C2 lines. The advantage of this feature is that the C2 
blends introduce a large asymmetry in the profile, allowing thus 
to readily constrain their contribution. Another advantage is that 
there are other nearby C2 lines that can be used to calibrate the 
blending C2 lines, and therefore to determine a reliable solar 
oxygen abundance. 
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Fig. 1. Observed solar spectrum (Delbouille et al. 1973) between 
5575 and 5580 A (red circles) and the best fit employing the 
<3D> model (solid blue line). The locations of the most promi- 
nent C2 lines used to estimate the C abundance are denoted with 
vertical lines, while the position of the [O i] line is marked with 
an arrow. The resulting solar C abundance is logec - 8.42 
with which the 5577.3 A feature implies a solar O abundance 
of log eo - 8.73. 



2. Analysis 

We used the Delbouille et al. (1973) disk-center solar spectrum 
(often called the Liege solar atlas), the Brault & Neckel (1987) 
disk-center intensity atlas (see Neckel 1999), and the Kurucz 
et al. (1984) solar- flux spectrum. These solar spectra are of ex- 
tremely high resolution (resolving power R ^ 500 000 - 700 000) 
and have a S/N of several thousands at 5577 A. The solar oxygen 
abundance obtained with these different atlases agree within cr 
= 0.04 dex. In Fig. 1 the Delbouille et al. (1973) solar spectrum 
around the 5577.3 A [Oi] line is shown as open circles. 

The oscillator strength for the [O i] 5577 A line is well deter- 
mined: the mean of the available calculations is log gf = -8.25 + 
0.03 (Galavis et al. 1997; Baluja & Zeippen 1988; Fischer & 
Saha 1983; Nicolaides & Sinanoglu 1973). We adopted log gf = 
-8.28 from Galavis et al. (1997), which is on the same scale as 
the log value adopted for the 6300, 6363 A forbidden lines 
(Asplund et al. 2004) within 0.007 dex; had we instead opted 
for the mean value of the above-mentioned theoretical values 
our derived O abundance would be 0.03 dex lower. The excita- 
tion potential of the [Oi] 5577 A line is 1.967 eV. According 
to the NIST database 3 the observed and Ritz wavelengths are 
5577.34 and 5577.339 A, respectively. Small errors in the wave- 
length translate into only small uncertainties (at the level of 0.01 
dex) in the O abundances. 

We included in the spectral synthesis other atomic and 
molecular (C2 Swan and CN red systems) lines present around 
5577 A, but the main contributors to the 5577.3 A feature are 
mostly three lines: the [Oi] line and the Pi27 and P226 lines 
of the C2 (1-2) band. Atomic lines were taken from the latest 
Kurucz line list^ and the CN line list described in Melendez 
& Barbuy (1999) was adopted. The line list for the 0-1 and 1-2 
bands of the C2 Swan system was constructed in a similar way to 
the 0-0 band described in Melendez & Cohen (2007), but with 
the line positions taken primarily from Tanabashi et al. (2007) 





<3D> model, log 



.605, log = 8.455 



' http://physics.nist.gov/PhysRefData/ASD/index.html 
^ http://kurucz.harvard.edu/ 
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Fig. 2. Upper panel: An enlargement of the observed (red cir- 
cles, Delbouille et al. 1973) and <3D>-based theoretical profiles 
(solid blue line) shown in Fig.[T] The dotted lines show the sepa- 
rate contribution of [O i] and other lines labeled accordingly. The 
5577.54 A P325 C2 line is blended with an unknown feature, as- 
sumed here to be Fe I. Lower panel: Same as above but assuming 
instead that the red feature is entirely due to P325 C2. Applying 
this 0.03 dex higher C abundance to the [O i]-i-C2 blend implies a 
smaller left-over contribution for [O i] and thus a 0. 1 3 dex lower 
O abundance at the expense of a significantly poorer profile fit. 



and complemented with Phillips & Davis (1968) for the higher 
excitation lines not included by Tanahashi et al. . Relative oscilla- 
tor strengths for the 0-1 and 1-2 bands were taken from Kokkin 
et al. (2007), normalizing those values to the laboratory oscil- 
lator strength of the 0-0 band recommended by Grevesse et al. 
(1991), /oo — 3.03 xlO"^. The rotational strengths were com- 
puted following Kovacs (1961) and the rotational dependence of 
the band strengths was taken into account (Dwivedi et al. 1978). 

A set of six different model atmospheres was employed: 
a three-dimensional hydrodynamical model of the solar atmo- 
sphere (Asplund et al. 2000; here denoted 3D) and its temporal 
and spatial average (here: <3D>) (Asplund et al. 2004), MARCS 
model (Asplund et al. 1997), Kurucz overshooting (Castelli et al. 
1997) and the latest no-overshooting (Castelli & Kurucz 2004) 
models and the semi-empirical Howeger-Miiller (1974) model. 
The spectral line formation was calculated in 3D and ID using 
the same spectrum synthesis code as described in Asplund et al. 
(2004). The ID calculations were also performed using the 2002 
version of MOOG (Sneden 1973); the different ID calculations 
with different codes agree at the level of 0.01 dex. 
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The contribution from the two C2 lines (Pi 27 and P226) to 
the 5577.3 A feature can be well constrained from the red asym- 
metry in the profile. In addition there are many neighboring C2 
lines of very similar excitation potential and line strengths and 
thus basically identical line formation. Note that the error in 
the relative strengths of the C2 lines is negligible, as the rota- 
tional strengths are precisely given by their Honl-London fac- 
tors (Kovacs 1961). Fig.[T]shows the best fit to the 5575-5580 A 
region of the Liege disk-center atlas based on the 3D-averaged 
model atmosphere; the locations of the most important C2 lines 
are marked with vertical lines. Clearly the overall agreement is 
quite satisfactory. Also encouraging is that the thus estimated 
solar C abundance is as expected for the different model atmo- 
spheres (Asplund et al. 2005). As seen in Fig. |2] the profile of 
the [Oi]-hC2 feature is very well described with this C abun- 
dance and a solar oxygen abundance of logeo = 8.7^ with 
the <3D> model. The resulting disk-center intensity equivalent 
width for only the [Oi] line is 1.23 mA, which corresponds to 
logeo = 8.72 for the full 3D model. The Kurucz et al. (1984) 
solar-flux and the Brault & Neckel (1987) intensity spectrum 
gave an O abundance 0.01 dex higher and -0.07 dex lower, re- 
spectively. Results for different model atmospheres are similar, 
as shown in TablelT] thus the [O i] 5577 A line is almost insensi- 
tive to the adopted model atmosphere. 

It should be noted that the predicted profile for the neighbor- 
ing feature at 5577.55 A is too weak with the C abundance from 
the overall fit to all C2 lines in the wavelength window 5575- 
5580 A (Fig.|2|i, suggesting that it is blended. Scouring available 
atomic databases did not reveal any likely candidates. We have 
therefore tentatively assigned the blending line as being due to 
Fe I but the exact choice is unimportant. Relying solely on the 
P325 transition and without accounting for this blend, the esti- 
mated C abundance would have been 0.03 dex higher Because 
of the larger contribution from C2 the [Oi] part of the 5577.3 A 
feature must decrease correspondingly, leading to a much re- 
duced solar O abundance: logeo ~ 8.6. In this case, however, 
the profile fit is very unsatisfactory, in particular in the red wing, 
as evident in Fig. |2] 

The main errors are the uncertainties in the log gf value of 
the [O i] line (assumed here to be 0.03 dex, which is the cr be- 
tween the different calculations; c.f. Galavis et al. 1997; Baluja 
& Zeippen 1988; Fischer & Saha 1983; Nicolaides & Sinanoglu 
1973) and the observational error. Part of the observational er- 
ror (0.04 dex) may be due to variation with solar activity (e.g. 
Livingston et al. 2007), and it was estimated by using differ- 
ent solar atlases. Note that Caffau et al. (2008) has also found a 
similar scatter (0.03 dex) between the [Oi] 6300A oxygen abun- 
dances obtained from different solar atlases. The other part of the 
observational error (0.05 dex) was estimated by employing dif- 
ferent assumptions with respect to the continuum determination 
and varying by 0.0 1 A the wavelengths of the C2 linefl and of the 
unknown feature at 5577.55 A. Given the excellent agreement 
between the oxygen abundance obtained from different model 
atmospheres for this spectral feature, we estimate that the er- 

^ Note that the abundance quoted here is only from the Liege inten- 
sity atlas (Figs. 1 and 2). The final abundance given in Table 1 is based 
on the mean from the three different solar atlases. 

According to Tanabashi et al. (2007) the line positions for those C2 
lines are not very precise, hence we allowed small changes by up to 
0.01 A in wavelength. Also, note that Allende Prieto & Garcia Lopez 
(1998) have shown that the Delbouille et al. solar disk-center intensity 
atlas have slight errors in its absolute wavelength calibration, although 
probably not over a wavelength scale of a few A, which is relevant here. 



rors introduced by problems in the model atmosphere should be 
0.04 dex at most. Non-LTE effects for this [Oi] line should be 
negligible (e.g. Altrock 1968; Takeda 1994). Adding the squared 
errors (0.03 dex in the log g/ value, 0.04 dex in the apparent line 
variation, 0.05 dex due to uncertainties in the continuum and line 
positions, and 0.04 dex in the models), we estimate a conserva- 
tive total error of 0.08 dex. 

3. Solar oxygen abundance from different 
diagnostics 

Asplund et al. (2004) have determined the solar oxygen abun- 
dance from the [O i] 6300, 6363 A forbidden lines, O i permitted 
lines, and the pure rotation and fundamental vibration-rotation 
OH lines, employing a 3D hydrodynamical simulation of the so- 
lar atmosphere as well as a ID MARCS (Asplund et al. 1997) 
and the Holweger & Miiller (1974) semi-empirical model atmo- 
spheres. Melendez (2004) has added oxygen abundances deter- 
mined from the first-overtone OH lines, as well as computed the 
solar oxygen abundance for the above features using the <3D> 
and Kurucz (Castelli et al. 1997) convective overshooting mod- 
els. Here we add the 5577.3 A forbidden line for the five model 
atmospheres described above, and we also have determined oxy- 
gen abundances for the latest Kurucz (Castelli & Kurucz 2004) 
model without convective overshooting. 

The results of the different computations are summarized in 
Table 1, where the mean oxygen abundance for the forbidden 
lines now includes the 5577.3 A line. The mean value given for 
the infrared OH lines includes all the OH mentioned above, but 
given half-weight to the pure rotation and the first-overtone line, 
because the pure rotation lines are quite sensitive to the detailed 
structure of the model atmosphere, while the first-overtone lines 
are very weak and the observational error is high. For the Oi 
lines we adopt the mean non-LTE results of Allende Prieto et 
al. (2004) and Asplund et al. (2004), with and without includ- 
ing inelastic H collisions, respectively. As discussed by Asplund 
(2005), the H collisions are based on the classical but highly 
uncertain formula by Drawin (1968), yet the O abundance with 
(and without) collisions is only k 0.03 higher (lower) than the 
mean O abundance adopted here. We also show in the last row 
the averaged solar oxygen abundance for each model atmo- 
sphere. In obtaining this mean value, half weight has been given 
to the result obtained from the molecular lines, because of their 
strong model atmosphere sensitivity. 

As can be seen, all model atmospheres favor a low solar oxy- 
gen abundance: logeo - 8.7. The reduction in the solar oxy- 
gen abundance from the historical value of ^ 8.9 (Anders & 
Grevesse 1989j3 to the present level is therefore not primarily 
driven by the use of a 3D hydrodynamic models. In this con- 
text, improvements in the input atomic and molecular data, more 
realistic non-LTE calculations and a more careful treatment of 
blends are equally important. The average O abundance for the 
3D model presented here is 8.67, only 0.01 higher than the value 
recommended by Asplund et al. (2004, logeo = 8.66 + 0.05). 
Note that the recommended O abundance obtained with the 
Holweger & Miiller (1974) model is only 0.04 dex higher than 
the 3D-based abundance but show the largest scatter (0.10 dex), 
mainly due to the large difference between OH and O i. 

' It should be noted that the solar O abundance of 8.83 given in the 
compilation of Grevesse & Sauval (1998) is only a preliminary analysis 
based on a Holweger-Miiller (1974) model with a modified temperature 
structure to remove various abundance trends with excitation potential 
and line strengths for species such as Fe i, OH and CH. 
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Table 1. Solar oxygen abundance from different model atmospheres. 



lines 


3D 


<3D> 


Holweger-MUller 


MARCS 


Kurucz overshooting 


Kurucz no overshooting 


[Oi] 5577 a;;. 


8.70±0.08 


8.71 


±0.08 


8.73: 


±0.08 


8.70: 


fc 0.08 


8.74 = 


t 0.08 


8.70 ± 0.08 


[Oi] 


8.69 + 0.05 


8.73 


+ 0.05 


8.75 


±0.05 


8.71 


±0.05 


8.77 = 


t 0.05 


8.73 ± 0.05 


Oi 


8.67 + 0.05 


8.68 


+ 0.05 


8.64 


±0.08 


8.72 


±0.05 


8.67 = 


t0.05 


8.64 ± 0.05 


OH 


8.61 +0.10 


8.68 


±0.10 


8.85 


±0.10 


8.73 


±0.11 


8.80 = 


t 0.10 


8.72 ±0.10 


Average 


8.67 + 0.04 


8.70 


±0.04 


8.71 


±0.10 


8.72 


±0.04 


8.73 = 


t0.07 


8.69 ± 0.05 



" Average abundance from the Delbouille et al. (0.02 higher), Brault & Neckel (0.04 lower) and Kurucz et al. (0.03 higher) atlases. 



Using 1/cr^ of the mean oxygen abundances given in Table 
[T| as weights (assuming a minimum cr - 0.04 dex) implies a 
weighted mean solar oxygen abundance from the six solar model 
atmospheres of log eo - 8.70 (cr = 0.02). 

Higher solar O abundances have recently been advocated by 
Caffau et al. (2008) and Ayres (2008) based on the atomic lines 
and an alternative 3D solar model computed with the C05B0LD 
code (Freytag et al. 2002); neither study includes the here em- 
ployed [Oi] 5577 A line. Caffau et al. (2008) finds log eo = 8.76 
with the main differences being the adopted equivalent widths 
and the less severe non-LTE abundance corrections for O i due 
to their inclusion of inelastic H collisions close to the Drawin 
(1968) recipe; it appears that the results of Asplund et al. (2004) 
and Caffau et al. (2008) in most cases are in excellent agreement 
when otherwise using the same line input data, in spite of some 
differences in the mean temperature stratification of the two 3D 
models. Ayres (2008) found log eo - 8.81 based on the [Oi] 
6300 A line, which is blended by a Ni i line (Allende Prieto et al. 
2001). In order to achieve the best overall profile fit, he had to 
use a solar Ni abundance lower by 0. 15 dex than the value given 
by Grevesse et al. (2007). This is well outside the quoted uncer- 
tainties in the laboratory gf -value (Johansson et al. 2003) and 
the solar Ni abundance. Clearly more work is needed in order 
to resolve the existing issues with solar O abundance determi- 
nations using the different indicators. We are currently working 
towards this goal using a new, further improved 3D solar model 
with a more realistic radiative transfer treatment (e.g. opacity 
sampling instead of opacity binning). 



4. Conclusions 

A low oxygen abundance is obtained from the [Oi] 5577.3 
A line, almost independent of the adopted model atmosphere: 
log eo = 8.71 ± 0.02 + 0.07 (0.02 dex is the cr from using dif- 
ferent models, and 0.07 dex is the error due to uncertainties in 
the the log gf value (0.03 dex), apparent line variation (0.04 
dex) and uncertainties in the continuum and line positions (0.05 
dex)). This value is close to the results of Asplund et al. (2004) 
and Melendez (2004) for other solar oxygen abundance indica- 
tors (Oi, [Oi] and OH lines). Including those abundances and 
employing six different 3D and ID model atmospheres, we esti- 
mate a solar O abundance of log eo - 8.70 (cr - 0.02) with only 
a small sensitivity to the employed model atmosphere. 

The low solar oxygen abundance that we propose here as 
well as the downward revision of other solar chemical abun- 
dances (Asplund 2000, 2005, 2007; Asplund et al. 2000, 2005, 
2006) pose a challenge to standard solar models in light of he- 
lioseismological observations (e.g. Antia & Basu 2005, 2006; 



Bahcall et al. 2004, 2005, 2006; Basu & Antia 2004, 2008; 
ChapHn et al. 2007; Guzik et al. 2005; Yang & Bi 2007). 
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